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Abstract 

We investigate the universality of microscopic eigenvalue correlations for Random Ma- 
trix Theories with the global symmetries of the QCD partition function. In this article 
we analyze the case of real valued chiral Random Matrix Theories [jS = 1) by relating 
the kernel of the correlations functions for (3 = 1 to the kernel of chiral Random Matrix 
Theories with complex matrix elements (/5 = 2), which is already known to be universal. 
Our proof is based on a novel asymptotic property of the skew-orthogonal polynomials: an 
integral over the corresponding wavefunctions oscillates about half its asymptotic value in 
the region of the bulk of the zeros. This results solves the puzzle that microscopic univer- 
sality persists in spite of contributions to the microscopic correlators from the region near 
the largest zero of the skew-orthogonal polynomials. Our analytical results are illustrated 
by the numerical construction of the skew-orthogonal polynomials for an probability 
potential. 
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1 Introduction 



Since its first application to the spacing of nuclear resonances Random Matrix The- 
ories (RMT) have been very successful in explaining the statistical properties of spec- 
tra. Originally, the so called Wigner-Dyson ensembles of matrices with independently 
distributed Gaussian matrix elements were introduced to replace the unknown nuclear 
Hamiltonian. More recently, the applicability of RMT has been related to the chaotic 
dynamics of the corresponding classical system and spectra of many chaotic systems with 
only a few degrees of freedom have been successfully described by RMT 0, |, ^. The 
successes of RMT has raised the question whether the statistical properties of spectra are 
universal. This question has been investigated in great detail within the context of RMT. 
The idea is to show that large deformations of the probability distribution of the random 
matrix elements leave the properly rescaled spectral correlations unaffected, whereas the 
average spectral density changes on a macroscopic scale. This program has been carried 
out most completely for the Hermitian random matrix ensembles (denoted by the Dyson 
index /5 = 2; for recent reviews see 0, ^) which are mathematically much simpler than 
real or quaternion-real random matrix ensembles (with Dyson index (3=1 and /3 = 4, 
respectively). Nevertheless, several universality proofs are available for these ensembles 
as well I, 0, 1, 1, [ig, ini, m, H, 0. 

In addition to the Wigner-Dyson ensembles there are seven other classes of Random 
Matrix Theories. They can be classified according to the Cartan classification of sym- 
metric spaces [|^. In this article we are interested in chiral Random Matrix Theories 
(chRMT). These are ensembles of random matrices with the chiral symmetry of the QCD 
Dirac operator jT^, The nonzero eigenvalues of these ensembles occur in pairs ±A. 
Therefore, A = is a special point, and the average spectral density on the scale of the 
average level spacing shows universal properties. With the average spacing of the eigen- 
values given by n/HN (with the total number of eigenvalues and S a parameter known 
as the chiral condensate), the microscopic spectral density is defined as [ITB 



Both Ps{u) and the microscopic fc-point correlation functions are universal. This has been 
shown in great detail for the chiral Unitary Ensemble (chUE), which is the ensemble of 
Hermitian chiral random matrices with no anti-unitary symmetries ||TB|, |T^, ^ ^, 
23], |2|, ^, g6|, |2§ [gl, ig, |T], |2|, g. The chiral ensembles with real or quaternion 



real matrix elements (known as the chiral Orthogonal Ensemble (chOE) and the chiral 
Symplectic Ensemble (chSE), respectively) are mathematically much more complicated. 
The general result for the microscopic spectral density of the chiral Gaussian Orthogonal 
Ensemble (chGOE) and the of the chiral Gaussian Symplectic Ensemble (chGSE) 
was first obtained by an explicit construction of the corresponding skew-orthogonal 
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polynomials ||3^, |3^ (several special cases were analyzed in ^). This method 

does not seem to be easy generalizable to non-Gaussian probability potentials. 

Recent progress was made by relating the kernel for the correlation functions of the 
chOE to the universal kernel of the chUE |TD[. This method was based on a generaliza- 



tion of an operator construction of skew-orthogonal polynomials for the Wigner-Dyson 
ensembles |0| to the chiral ensembles. Remarkably, the skew-orthogonal polynomials do 
not enter in the relation between the kernels. Indeed, an elegant construction of the cor- 
relation functions relying only on operator relations was recently given in [|10|, |ll]]. The 
operator method was successfully applied to Gaussian Orthogonal and Gaussian Sym- 
plectic Ensembles with an additional fermion determinant [Q. Universality for the so 
called massive chiral ensembles for j3 = 1 and /3 = 4 was shown by relating them to 
the corresponding massless ensembles [|14|. In this way analytical results for the massive 



spectral correlators could be obtained |14, 41 



The universality proof based on the relations between kernels given in |JTO| is incomplete 
for (3 = 1. The reason is that the correlation functions for /3 = 1 depend on a integral over 
the positive real axis of a derivative of the kernel for (3 = 2. Naively, one would expect 
that the presence of such nonlocal contributions would lead to nonuniversal behavior. The 
resolution of this puzzle, and thus the completion of the universality proof for /5 = 1, is 
the primary objective of this paper. As was already discussed in for /? = 4 there is 
no such problem. As secondary objective, we illustrate some of the arguments given in 
T0|] by the analysis of the asymptotic behavior of the skew-orthogonal polynomials for a 



quartic probability potential. 

The results of this article are relevant for Dirac spectra of QCD with two colors in 
the fundamental representations (or for Dirac spectra of QCD with staggered fermions 
in the adjoint representation). Indeed, the microscopic spectral density of the chGOE 
has been observed in lattice QCD ^ and in instanton liquid simulations |Q. Many 
more results justifying the chiral Random Matrix description of the microscopic spectral 
density in lattice QCD have been obtained for /3 = 2 and /5 = 4 [HH HBI, H^, iS, EUf (for 



a recent review and a complete list of references see Our results also apply to the 

superconducting ensembles with Dyson index (3 = 1 (with joint eigenvalue density given 
by a special case of the chiral ensembles) as well as to two-sub lattice models [0, . 

This paper is organized as follows. Chiral Random Matrix Theory is introduced in 
section 2. In section 3 we discuss the relation between the kernels for (3 = 1 and (3 = 2. 
Most of this section already appeared in |T^. Two examples, the Gaussian case and the 
quartic probability potential are worked out in detail in section 4. In section 5 we derive 
a novel asymptotic property of the skew-orthogonal polynomials, allowing us to complete 
the proof of This is the most important result of this paper. Universality of the 
microscopic correlations is shown in section 6. In this section we also give the explicit 



3 



universal expressions for the microscopic spectral density and the microscopic kernel which 
allows to calculate all correlation functions. In section 7 the skew-orthogonal polynomials 
for a quartic probability potential are studied numerically and concluding remarks are 
made in section 8. 



2 Chiral Random Matrix Theory 

In this section we introduce a Random Matrix Theory with the global symmetries of 
the QCD partition function. We first define the Dyson index /? of a Dirac operator. Its 
value is equal to the number independent degrees of freedom per matrix element and 
is determined by the anti-unitary symmetries of the Dirac operator. Any anti-unitary 
symmetry operator can be written as f/ = AK with A unitary and K the complex 
conjugation operator. The operator f/^ is unitary and is thus proportional to the identity 
in an irreducible subspace of the unitary symmetries of the Dirac operator. One can easily 
convince one- self that the only two possibilities are A? = \ and A? = —\m. this subspace. 
These two cases are denoted by the Dyson index /3 = 1 and /5 = 4, respectively. If there 
are no anti- unitary symmetries, the Dyson index of the Dirac operator is /5 = 2. For 
/3 = 1 it is possible to find a basis for which the Dirac matrix is real for all gauge field 
configurations. For /5 = 4 it is possible to construct a basis for which the Dirac matrix is 
quaternion real for all gauge field configurations. For /3 = 2 the matrix elements of the 
Dirac operator do not have any reality properties and are given by complex numbers. The 
Dyson index of the Dirac operator depends on the representation of the gauge fields and 



may be different for the discretized and continuum versions of the Dirac operator |52 
For example, the continuum Dirac operator for Nc = 2 and fermions in the fundamental 
representation is in the class /? = 1, whereas the corresponding staggered Dirac operator 
is in the class /3 = 4. For gauge fields in the adjoint representation the Dyson index of 
the continuum Dirac operator is /3 = 4, whereas the Dyson index of the staggered lattice 
discretization is /5 = 1. 

In a chiral basis in the sector of topological charge v the Dirac matrix has the block 
structure 

/ n ri \ 

(2) 

where C is an n x (n + z/) matrix with reality properties given by the Dyson index. For 
generic values of its matrix elements D has exactly v zero eigenvalues. In QCD, the matrix 
elements of D depend in a complicated way on the gauge fields, which are distributed 
according to the QCD partition function. In chiral Random Matrix Theory (chRMT) 
we replace the matrix elements of the Dirac operator by space-time independent random 
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numbers with probability distribution given by the partition function 

Z{m) =m''J dCdet^'fiD + m)e-^T^^(^'^). (3) 

Here, Nf is the number of quark flavors with mass m (below we only consider the mass- 
less case 171 = 0). The integral is over the independent degrees of freedom of the matrix 
elements of C. The probability potential V{x) is in general an arbitrary polynomial 
V{x) = OfcX^ (with Op > 0). For the Gaussian chiral ensembles, which are mathe- 
matically much simpler, the potential is V{x) = S^x. In that case the parameter S is 
related to the average spectral density 

P(A) = (E<5(A-A,)) (4) 

k=l 



by the Banks-Casher formula [53 



E=n.n^42^. (5) 

The prime indicates that the argument of p(A) should be near zero but much larger than 
the smallest eigenvalue. For this reason S is interpreted as the chiral condensate, the order 
parameter of the chiral phase transition. The integrals in (|^) can be easily evaluated in the 
thermodynamic limit. The result coincides with so called finite volume partition functions 



which were first derived on the basis of chiral symmetry |55| . 

In this paper we will study the chiral Orthogonal Ensemble (chOE) and show that 
the microscopic spectral density does not depend on the coefficients of the probability 
potential. Our starting point is the joint probability distribution of the eigenvalues of 
the Dirac operator. It can be obtained form the partition function by making a polar 
decomposition of C 

C = UAV-\ (6) 

with U and V orthogonal, unitary or symplectic matrices for (3 = 1, f3 = 2 and /3 = 4, 
respectively, and A a semi-positive definite diagonal matrix. For the massless case the 
joint probability distribution in terms of the squares of the eigenvalues, Xk = A|, is given 
by 

p(.„...,x„) = |A({x.})|^n-^"'"''^''"''V^i:.^(-) (7) 

k 

where the Vandermonde determinant is defined by 

A{{x,}) = Y[{xk-xi). (8) 

k<l 
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The exponent of the Xk for /3 = 1 will be denoted by 

a = Nf -l + \v\/2 + l/2, (9) 
and below we consider the joint probability density 

p(xi,---,xO = |A({x,})|'^ne-^^"^''=^- (10) 



k 



with probability potential given by 



—V{x) — alogx. (11) 



For technical reason we will restrict ourselves to even n and use the notation n = 2n. 

The spectral density, and in general the /c-point correlation functions, can be obtained 
from the joint eigenvalue density by integrating over all but k eigenvalues. For (3 = 2 this 
can be simply done by exploiting the orthogonality of the polynomials 



/ dxe-^^-^''^Pl\x)Pf''{x) = Ski. (12) 
Jo 

The resulting spectral correlation functions can be expressed in terms of the kernelQ 

n-l 

K'rr{x,y) = T.Pti^)Ptiy)- (13) 

k=0 

Universality can be established by showing that in the microscopic limit the so called 
wave functions P|"(x) exp(— 0a(x)) depend only on the potential through a scale factor 



determined by the average spectral density |19|. With E = 7rp(0)/2n the microscopic 
limit of the kernel for the chUE is given by (a factor 2^Jm} from the integration measure 
has been included) 

lim „o,.o A„ — = Juv 



(14) 

This kernel is known as the Bessel kernel |5^. Sometimes it is simpler to use an integral 



representation of the Bessel kernel given by 

B^''{U,V) = f\dtJ2a{ut)J2a{vt). (15) 

Jo 



^For later convenience we do not follow the usual convention to include the weight functions in the 
kernel. 
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3 Relation between the kernel for P = 1 and (3 = 2 



For the orthogonal ensembles the integrations over the eigenvalues can be performed by 
means of orthogonality relations for the skew-orthogonal polynomials of the second kind 
1 3^, 0. These polynomials are defined by the scalar products 



{Rk-, Ri)r — Jkh (16) 

with the nonzero matrix elements of Jki given by J2k,2k+i = —J2k+i,2k = —1- The skew- 
scalar product is defined by 



{f,g)R= / dxe-'*^^^^f{x)Zg{x), (17) 
Jo 

where we have introduced the operator Z by 0, [10|, 



Zg{x)= dye'^^^''h{x -y)e-^^^y'^g{y). (18) 
Jo 

As usual, e{x) = x/2\x\. All correlation functions can be expressed in terms of the kernel 







with the pre-kernel defined by 



2n-l 

kl{y,z)= Y: Rt{y)J^,R';{z). (20) 

i,j=0 

For example, the spectral density is given by 

p{x) = KI{x,x)-\k^^{^,x). (21) 

We construct the skew-orthogonal polynomials for the chUE by a generalization of 
an operator method ||T0| introduced by Brezin and Neuberger The skew-orthogonal 
polynomials are expressed in terms of orthogonal polynomials p^a+fc ^j^j-^ weight function 
x2«+^exp(-ny(x)/2), 

R'i{^) = iZT,3Pr^\x). (22) 

j=0 

We will derive recursion relations for the expansion coefficients Tij. It is useful to introduce 
the operators 

X, X^d, X'^Z, (23) 
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and the operator 

L = X^d- Xy(X) + (24) 

The coordinate operator and the derivative operator are defined by Xf{x) = xf{x) and 
df{x) = f'{x), respectively, and the operator L is the inverse of X~^Z, i.e., 

LX''Zg{x)=g{x). (25) 

The matrix elements of the operators in (^) and ( ^4]) with respect to the basis P^°'^^{x) 
will be denoted by Xm, D^i, Ym and L^i, in this order. For an operator A the matrix 
elements are defined by 

Aki = / rfa;e-2<^"(")+^^°*5"Pf"+^(x)APf +^(x). (26) 
Jo 

By partial integration it can be shown that L and D are related by 

Lki = ^iDM-Dik). (27) 

For integer values ofQ b > 1 the matrix elements of L thus vanish for |fc — Z| > max(6 — 
1, 6 + p — 1) (where p is the order of the polynomial probability potential). The optimum 
value of b is thus 6=1 which will be our choice in the remainder of this article. 
The orthogonality relation (|TB|) can be written as 



TYT^ = -J. (28) 
By acting with L and T^JT on YT^ it follows that 

L = T^JT. (29) 
In matrix form this equation can be rewritten as 

oo 

Lkl = y^[^2p+l.fc^2pj — r2p,fc^2p+l,/]- (30) 
p=0 

For known L^i the coefficients Tj^ can be determined recursively from this relation. If the 
vanish outside a band \k — l\ > p (as is the case for V{x) given by a polynomial of 

order p), the coefficients Tik are nonzero only inside a band i — k < M < 2p. 
The pre-kernel can be written as 

2n-l 

i,j=0 k<i l<j 
2n-l 

= T.T.Pk''^\^)TlJ^,T,iP^''^\y)-R{x,y), (31) 

k,l=0 i,j 



^For 2a — the matrix L remains a band matrix even for 6 = 0. 
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where remainder term is given by 



2ra+Af-2 
i,j=2n k<i l<j 



(32) 



There are now no restrictions on the summation over i and j in the last hne of (|3l|) and 



the relation (pOf ) can be used to simplify the expressions. For finite M and a smooth 
dependence of the coefficients Tij on the order of the skew-orthogonal polynomials the 
number of terms in R{x,y), and thus the contribution of R{x,y) to the pre- kernel, is 
subleading in 1/n. In the next section, this will be shown explicitly both for a Gaussian 
and a quartic probability potential. To leading order in 1/n we thus find 



2n-l 



k,l=0 
1 2n-l 



(33) 



k,l=0 



By re-expressing the matrix elements of D in terms of the operators xd^ and ydy we 
find the following relation between the kernel for the chiral Orthogonal Ensemble and the 



kernel for the chiral Unitary Ensemble |]10 



kRix,y) = ^{ydy - xd,j,)Kl^^^{x,y). 



(34) 



This relation was first obtained for the Gaussian case in ||3^. Since it has been shown 
that K'^^^{x,y) is universal [19|, we thus have proved that the pre-kernel is universal 
The only problem is that an integral of the pre-kernel over the complete spectrum 
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contributes to the spectral density and the spectral correlators. Even in the microscopic 
limit this results in contributions from non-universal regions. Before going to the general 
case, we first analyze in detail the chGOE and the case of a quartic probability potential. 



4 Two Examples 

In this section we study the quadratic and the quartic probability potential. In the 
first case, the matrix elements of L and the skew-orthogonal polynomials will be derived 
exactly, whereas in the second case only asymptotic results for large order polynomials 
will be obtained. 



4.1 The chiral Gaussian Orthogonal Ensemble 

In this section we study the chiral Gaussian Orthogonal Ensemble by means of the op- 
erator construction discussed in the previous section. The weight function is given by 
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w{x) = x^^^^e and the corresponding orthonormal polynomials are the Laguerre 
polynomials, 

with normalization constants 

and ft;: = £(^±fl±i). (36) 



" n"+i A;! 
The matrix elements of L follow immediately from the recursion relation 

xd^Ll^^\x) = kLf^\x) -{k + 2a + l)Ll''_\\x), (37) 

and are given by 



Lki = i[(/ + 2a + l) 



TSi,k-i]. (38) 



One easily verifies that the recursion relation (^) does not have a solution for diagonal 
matrices T^. A solution is obtained by taking T2k,2k, T2k+i,2k+i, T2k+i,2k and 72fc+i,2fc-i 
as the only nonzero coefficients. In the normalization R2k{x) = (nx) / \J (i.e. 
T2k,2k = 1) the recursion relations (pOf) simply read 



^2fc+l,2fc + l — -^2fc+l,2fc; r2fe+l,2fe-l — -^2A;,2fc-l; (39) 



and the skew-orthogonal polynomials are thus given by 



-R2fc(^) — / „ , , LlV'^i^^ 

'2A; 



(x) - + 2a + 2) ^ . (2A: + 1 + 2a) 1 

2 S2,+i 2 ^52a+l 

+ T2k+l,2kLlt+\x). (40) 

The coefficients T2k+i,2k are not fixed by the orthogonality relations. Indeed, this is the 
well-known property that the odd order skew-orthogonal polynomials are only determined 
up to a multiple of the even order polynomials of one degree lower. One can verify 
that these polynomials are normahzed according to (-^2^+15-^2^'''^) — 1; ^^d, with an 
adjustment of the normalization, they coincide with the polynomials obtained in [Rl] for 
a specific choice of coefficient T2k+i,2k- 

In fact, we can calculate the pre- kernel directly from (|3TD using the matrix elements 
of T^JT without relying on explicit expressions for the skew-orthogonal polynomials. 
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Because only one term contributes to R2k{x), there are no terms in (|30[) that are outside 
the summation range in (^l]). We thus have 

2n-l 1 -1 

kRix,yr= E Lf^\nx)L,iLr^\ny) = -{ydy - xd,)Kl^^\x,y), (41) 



which was first obtained in from the exphcit properties of the skew-orthogonal poly- 
nomials. Using recursion relations for the Laguerre polynomials the pre-kernel can be 
rewritten as 

-{ydy-xd.)K:{x,y) = -^(LU{nx)Lt\ny) - LU{ny)Lt-\nx)) 

^ fc=0 ^"^fc 

= l(9,-a,.)i^ri(x,2/), (42) 

where the factor k + a has been absorbed in the normalization of the orthogonal polyno- 
mials. 

According to (pT]), the spectral density is given by 

1 fX 1 roc 

- / rf2e-<^"(^)-'^-(^)(9, - 9^)i^f (x, z)-- dze-^^^'^-^'^^''\d, - d^)Kl\x, z) (43) 

2 JO 4 JO 

The microscopic limit of the first term results in an integral over the universal Bessel 
kernel. However, it is not possible to the interchange the integral and the microscopic 
limit in the second term. To see this we return to the definition of the pre-kernel. Then 
the second term is given by 

1 foo 2n—l 

- - / d2e"^"(") y R1{x)Ji,m{z)e-^'^^^\ (44) 
^ ■'^ i,3=0 

We thus consider the integral 

/•OO /"OO 

/ dze-'^'^^'^WAz) = / d;zz°e-""/2^f(z). (45) 

JO JO 

By using the explicit expressions for the skew-orthogonal polynomials and the relation 

n 

^r'(2x) = E Ll-m{x)Ll.{x), (46) 

m=0 

it follows that the integral over the odd order skew-orthogonal polynomials vanishes. The 
integral over the even order skew-orthogonal polynomials is up to a normalization constant 
given by 



[n 



/2)a+l 



j dzz''e-'''''^Lll+\nz) = 1. (47) 
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Let us now calculate the integral by interchanging the integration and the microscopic 
limit. Using the asymptotic relation relation between Laguerre polynomials and Bessel 
functions 

L^(nx) ~ a;-"/V«(2r2v^) (48) 
the microscopic limit of the integral is given by 

/ lim dzz'^e-'"'/^^^ Lt+\nz) = / dwJ2a+ii2w) = -. (49) 

Jo p ^ °° Jo 2 

p/n — fixed ^ 

Exactly half of the integral is missing. In the next section we will argue that this is a 
general feature of the skew-orthogonal polynomials. However, the integral over the odd- 
order polynomials does not vanish in general. The result for the microscopic spectral 
density is thus given by 

1 1 r°° 

Ps{u) = - / dwid^ - du)B^\u,w) - - / dw{d^ - du)B^\u,w), (50) 

2 JO 2 Jo 

where an additional factor of 2 has been included in the second term and the Bessel kernel 
is defined in (p!^). 



4.2 Chiral Orthogonal Ensemble with Quartic Potential 

In order to construct the skew-orthogonal polynomials using the Brezin-Neuberger for- 
malism 1^, we need an expression for the derivative of orthogonal polynomials. To derive 
such relation we start from the recursion relation 



xPk{x) = -rk{Pk+i - Pk) + Sk{Pk - Pk-i): (51) 
The coefficients and Sk are related by 

= (52) 

and hk = Jl^ dxw{x)P^{x) is the normalization integral. The recursion relation (^) is 
valid for orthogonal polynomials normalized according to Pk{0) = 1. To make contact 
with the analysis of we will use this normalization in this section. 

The derivative of the polynomials for arbitrary weight function is given by |^ 

^ Piix)P iiy) 

1=0 



yPkiy) = r dxw{x)j:^i^^xPUx) 

Jo t^r, hi 



kPk{y) + r dxw{x) V ^^^^4^^xP^(x) 
Jo r^r. hi 



.-0 

Pi{x)Pi{y) 



kPkiy) - r dxw'ix) Y: ^^^xPk{x). (53) 
Jo hi 
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where the terms following the last equal sign have been obtained by partial integration. 

Next we derive the asymptotic form of large order skew-orthogonal polynomials for 
a quartic probability potential. In terms of the Xk = Xl the weight function is given by 
w{x) = a;2a+ig-nx2/2^ Only the terms I = k — 1 and I = k — 2 are nonvanishing in the last 



sum in (|53D. Using the recursion relation ([51| ) to calculate the integrals we find 



yP'k{y) = kPk{y) - nsk[rk + Sk + Vk-i + Sk-i]Pk-i{y) + nskSk^iPk-2{y)- (54) 

Taking into account the normalization of the orthogonal polynomials the matrix elements 
of L are given by 



^kl 



1 1 

- / dxw{x)^==[Pi{x)xPl{x) - Pk{x)xP'i{x)] 
/ "'0 y/hkhi 

n jhk-2 . n Ihk-i r , , , n 

(>l,k-2^\l -^SkSk-l - Oi,fe-l-y -^Sk[rk + Sk + Tk^l + Sfc_iJ 



(55) 



As expected, the Lki vanish for |A; — / | > 2. 

In the limit n — oo the leading order contributions to the kernel are for terms with 
large values of k and /. The large-fc asymptotic behavior of the coefficients in the recursion 
relation ( ^T] ) was obtained in [n5| for an arbitrary polynomial probability potential (but 



for integer values of 2a). For fixed t = k/n, the continuum limit of the coefficients can be 
parameterized as 



For k — s> oo, only Sk and enter in the matrix elements of L. They are given by |T9 

r(t) = s{t) = y|. (57) 

The leading order asymptotic result for the matrix elements of L thus reads 

_ x 4nt Ant nt 

J^kl — Ol k-2— Oi k-l— h Ok,l-l— Ok^l-2 — - (,58) 

D D D D 

By inspection of the recursion relation (^) one easily finds that for a quartic potential 
(i.e. p = 2) the skew-orthogonal polynomials can be expressed as 

R2k{x) = P2k{x) +T2k,2k-lP2k-l{x), (59) 
R2k+l{x) = T2k+l,2k+lP2k+l{x) + T2k+l,2kP2k{x) + T2k+l,2k-lP2k-l{x) + T2k+l,2k-2P2k-2{x) ■, 

(60) 
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where the orthonormal polynomials have been denoted by Pk{x) = Pk{x)/^/hk- The 
coefficients Tj ^ obey the recursion relations 

^2fc-l,2fe-l ~ ^2fe,2A:-1^2fc+l,2fc-2 = -^2fc-l,2fc-2 5 ^2fc+l,2A:-2 = " -^2A;,2A:-2) 

^2A:,2/c-1^2/c+l,2fc ~ ^2fc+l,2fc-l = -^2A:,2/c-l; ^2fc,2A;-1^2A:+l,2A:+l = -^2fc+l,2fc-l- (61) 

Using the asymptotic values for the matrix elements of L we find a recursive equation for 
the coefficients T2p^2p-i, 

^2fc,2A:-1^2fc-2,2fc-3 + 4T2fc_2,2fc-3 + 1=0. (62) 

This recursion relation has two fixed points given by the roots of 

+ 4x + 1 = 0. (63) 

For large k the coefficients T2k,2k-i should depend smoothly on k and are thus given by one 
of the fixed points. The polynomial corresponding to the stable fixed point, x = —2 — -\/3, 
given by P2k{x) — (2 + V^)P2k-i{x), is negative for x = (we work in the convention 
Pi{0) = 1, and -P2A:(0) and -P2fc-i(0) are equal to leading order in 1/k) and positive for 
X — oo and thus has an odd number of zeros. Since even orthogonal polynomials should 
have an even number of zeros, the relevant solution is thus given by the unstable fixed 
point X = —2 + v^- This counter- intuitive result is a refiection of the numerical instability 
of the iterative construction of orthogonal polynomials. A numerical confirmation will be 
given in section 7. All other coefficients simply follow from the relations ( |6l|) resulting in 
the polynomials 

R2k{x) = P2k{x)-{2-V3)P2k-lix), 

I \ ""^^ n / \ 4nt ~ , , nt ~ , , 

R2k+Ax) = ■^^^—-^P2k+l{x) + —P2k-l{x) - —P2k^2{x) 

+ T2fc+i,2fc(P2fc(x)-(2- V3)P2fc-i(x)). (64) 

However, we do not need the explicit expressions for the skew-orthogonal polynomials. 
The pre-kernel can be expressed directly in the matrix elements of L and two additional 
terms that are outside the summation range in (pl|), 

2n-l ^ 

kR{x,y) = J2 rT-^ Pk{x)LkiPi{y) 
k,i=o V^kni 

^2n,2n-1^2n+l,2n-2 i 



+ T [-P2.~i(^)A.-2(y) - P2n-i{y)P2n-2{x) 

V il2n~lil2n-2 



(65) 



Because lim„^oo T2n,2n-i = — 2 + -\/3, the additional terms are of the same order of mag- 
nitude as each of the terms in the sum. Therefore, to leading order in l/n the relation 
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between the pre-kernel for (3=1 and the kernel for /3 = 2 is the same as for the Gaussian 
case. However, it is not clear whether the interchange of the integral and the microscopic 
limit in Kr{oo,x) also proceeds in the same way. This question will be analyzed in the 
next section for an arbitrary probability potential. 



5 A property of large order skew-orthogonal polyno- 
mials 

In this section we will derive an asymptotic relation for large order skew-orthogonal poly- 
nomials. Our starting point is the skew-orthogonality relation 

rco roc 

{Rk, Ro)r = dx dyRk{x)e-^''^^^-'^'^^yh{x - y) = 0. (66) 
Jo Jo 

For asymptotically large k we can distinguish three different domains in the integration 
over X. The region near a; = 0, the region around the largest zero of Rk{x), and the 
oscillatory region which is in between these two regions. We split the integration over 
X into two parts separated by M chosen to be inside the oscillatory region. In our 
normalization the spacing of the smallest eigenvalues scales as Ax ~ The main 

contribution to the integral over [0, M] is from the region around the smallest eigenvalues 
X ~ l/n"^, whereas the ^/-integral has contributions up to nV{y) ~ 1. If the potential 
behaves as with p > 1/2 near ?/ = 0, to leading order in l/n the contribution to the 
I/— integral is from the region with y > x and e{x — y) = —1/2. The main contribution 
to the integral over [M, oo) is from the region near the largest zero of Rk{x). To leading 
order in 1/n we can then replace e{x — |/) — > 1/2. The integrals over x and y factorize 
and we obtain the following asymptotic relation 



/ dxRi{x)e-'t"^^''^dx = / dxRl{x)e-^^^''^dx. (67) 
Jo Jm 

I dxRl{x)e-'*'^^''^dx = 2 / M^(x)e-^'"(^)rfa;, (68) 
Jo 



This implies that 



which is valid for A; ^ oo provided that the r.h.s. is independent of M in the oscillatory 
region. 

In our derivation we have made the assumption that the skew-orthogonal polynomi- 
als show the same oscillatory behavior as the regular orthogonal polynomials. This is 
certainly true if they can be expressed in a finite number of regular orthogonal polyno- 



mials which is the case for a finite order polynomial probability potential [|I0| . Of course, 
a different oscillatory behavior is possible if the leading order asymptotic terms in the 
skew-orthogonal polynomials cancel. This results in contributions that are subleading 
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in l/n. A priori, it cannot be excluded that the skew-orthogonal polynomials near zero 
and the l.h.s. of (|68| ) are of the same order in 1/n and thus both contribute to the ker- 
nel. For example, this is the case for the odd-order skew-Laguerre polynomials (0) with 
T2k+i,2k = 0. However, in this case the l.h.s. of (^) vanishes, and we do not have to 
worry about the asymptotic behavior of the r.h.s. of (0). In general, there is no reason 
to expect that the leading order asymptotic expansion of the P^""^^ cancels. For example, 
the asymptotic behavior of P^'^^^ in the oscillatory region depends smoothly on k, and 
from the explicit expressions for the even order skew-orthogonal polynomials in (^) it 
then follows that the leading order asymptotic behavior does not cancel. For the odd 
order skew-orthogonal polynomials a cancellation can only be achieved by fine tuning the 
coefficient T2k+i,2k- 

For a finite order probability potential, the generic asymptotic behavior of the R%{x) 
is thus the same as that of e~'^°-'^^^ P'^"'^'^ {x) . In the region near zero and in the oscillatory 
region, it is given by J2a+i{cVkx) / ^/x (with c a constant that can be obtained from 
the recursion relations). If the integral (^) is nonvanishing to leading order in 1/n, we 
thus find from the asymptotic behavior of e~'^''^^'' P^""^^ [x) that the integral converges for 
k —>■ oo and M inside the oscillatory region. 

As will be shown below, the even order skew-orthogonal polynomials are determined 
up to a multiplicative constant. From the orthogonality relations it is clear that the odd 
order polynomials R2k+i{x) are only determined up to the addition of a multiple of R2k{x). 
We may use this freedom to chose a normalization such that 

POO 

c/xe-'^"(^)i?^fc+i(x) = 0. (69) 

In this way we avoid the ambiguity in the asymptotic behavior of the odd order skew- 
orthogonal polynomials. 

The main contributions to the integral in the r.h.s. of (|68|) are from the region close 
to X = and the integrand can be replaced by its microscopic limit. Our final result is 

/•oo /-oo 

lim / i?^fc(x)e-'^"(^)rfx = 2 / lim Rlj^{x)e-^^^^^ dx. 



k/n — consz. xn^ — 2" 

k / n — const. 



(70) 



Interchanging the microscopic limit and the integral gives rise to an extra factor two. If 
the integrand in the kernel K'j;^{oo,x) is replaced by its microscopic limit, the same extra 
factor two has to be included, 

lim KUoo,x) = 2r lim du;e-^"("'/"')-*"(^/"')4't^(^, A)- (71) 



xn — 2 xn — z 



We emphasize that this relation is based on the asymptotic properties of the even order 
skew-orthogonal polynomials only. 
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6 Universality of the Chiral Orthogonal Ensemble 



The proof of universality of microscopic spectral correlation functions for the chiral or- 
thogonal ensemble for a finite polynomial probability potential is now straightforward. It 
is an immediate consequence of the following three results, i) The relation between the 
pre-kernel for the chOE and the chUE kernel is independent of the probability potential 
to leading order in 1/n [|10] . ii) The microscopic limit of the chUE kernel is universal if 
the eigenvalues are expressed in units of the average level spacing |jl9[. iii) Because of a 
novel asymptotic property of large order skew-orthogonal polynomials, the microscopic 
limit and the integrals that occur in the spectral correlation functions can be interchanged 
at the expense of a factor of 2. 

All fc-point correlation functions of the chOE can be expressed compactly as quaternion 
determinants of quaternions 

o/^ K^R,n{xk, Xl) - \K%Soo, Xl) X.? dz[K%,^{xk, z) - \K%Soo, z)] - e(x - y) ~ 



where K'^{x,y) is defined in ([191) . We already noticed that the spectral density is given 
by p{x) = Kf^^^XjX) — ^K'^^{oo,x). The two-point cluster function is given by @ 

Tix,y) = ^TTS{x,y)S{y,x). (73) 

The universal microscopic kernel for the chOE is obtained by taking the microscopic limit 
of the integrands and replacing the factors 1/2 — » 1, 

b[u,v) = hm ———Tb[ 



Q^iu, v) - Q%oo, v) T?N^ jS dw^lQ^^iu, w) - g"(oo, w)] 
^f^d^2Q''{u,v) (5"(t;,M) - Q"(oo,m) 



2\/uv 



(74) 



where Q"'{u,v) is the microscopic limit of the kernel K'^^{x,y). Universality then follows 
from the relation between the microscopic kernels for (3 = 1 and (3 = 2 and the universality 
of the microscopic limit of the kernel for (3 = 2. The universal microscopic result for 
Q'^iu^v) is given by 

1 Vj^ iP' 

= -/ d[w^]{wv)^''{d^2-d,2){wv)'^''''/^B^''{v,w). (75) 
4 JO 
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By using the integral representation ([T5|) of the Bessel kernel we find the following explicit 
representation of the kernel 



Q''(u,v) = - dw t^dt 
4 JO Jo 



uw 



J2a{uwt)J2a+l{vt) - J2a+l{uwt) J2a{vt) 



(76) 



One of the integrals can be performed analytically. By using identities for Bessel functions 
one ultimately finds the result 



2v 

J2a+l{v) 

2v 



dwJ2a+l{w) - 1 
dwJ2a-l{w) - 1 



(77) 



The microscopic spectral density for the chOE is given by Psiu) = 2u{Q°'{u, u)—Q°'{oo, u)). 



For its universal form we thus obtain 



Ps[U) 



y / t^dt 



dw [wJ2a{uwt)J2a+l{ut) - J2a+l{uwt) J2a{ut)] + -J2a+l{u) . 



(7J 



This expression for the microscopic spectral density can be simplified to 



Ps[U) 



u 

2 

u 

2 



J2a+liu) - J2a+2{u)J2a{u) 



- -J2a+l{u) 
Jlai.'^) - -ha+liu) J2a-l{u)\ - ]-J2a+l{u) 



dwJ2a+l{w) - 1 
dwJ2a-l{w) - 1 



(79) 



The simplified results for Ps{u) and the kernel were first obtained in [58| and was derived 
independently in 



(with a typo, see also [0, Q). The first term can be recognized 
as the microscopic spectral density for the chGUE |61|, The microscopic limit of the 
two-point correlation function follows immediately from (|73[) and the microscopic limit of 
the kernel. 

Recently, universal results for massive spectral correlators at /? = 1 and (3 = 4 have 
been obtained by relating the kernels for the massive correlators to the corresponding 
massless kernel . Similar relations have been derived for the Gaussian case ET] . 



7 Numerical study of the x -potential 

In this section we explicitly construct the skew-orthogonal polynomials for an potential 
and test the asymptotic results obtained in previous sections. We consider the distribution 
of the squared eigenvalues = A| on [0, oo( with weight function (pa{x) = — a log a; . 
For reasons of numerical accuracy we only consider integer values of a. 

The skew-orthogonal polynomials can be expanded in terms of monomials x'^ as 

Rtix) = rl,x' + rl,^,x'~' + ■■■ + rl,. (80) 
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Figure 1: Ratios of the integral /^(M) = /o^^/"''^'^^^' dxe-<^<'(^)i?fc(x) versus M for a 
Gaussian (left) and a quartic (right) probability potential. The parameter a and the 
value of k are given in the label of the figures. 
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The coefficients are determined from tlie ortliogonality relations (|T6D. They can be reduced 
to a system of hnear equations. For the polynomials of even order, R^^k 5 finds 

''00'2fc,0 "T ''0l'2fc,l -r ■ ■ ■ -|- ''0,2A:' 2fc,2fc ~ 

''10'2fc,0 "T ''ll'2fc,l -r ■ ■ ■ -t- ''l,2fc'2fc,2fc "'-'5 

''20'2fc,0 "T ''2l'2fc,l "T "T ^2, 2k' 2k, 2k ~ l^i-'-J 

^2A;-1,0' 2fc,0 "T ''2fc-l,l' 2fc,l "T " " " "I" ''2/c-l,2/c' 2fc,2/c ~ 

For the polynomials of odd order 2/c + 1, the first 2k — 1 equations have one more 
term, ^"2fc+i'"2fc+i,2fc+i5 the r2fc ,j are replaced by rg^+i j. The normalization equation 
{R2k+i,R2k)R = 1 reads 

'"2A:,2fc[^2fc,0^2A:+l,0 + '^2A:,1^2fc+l,l + ' ' ' + ^2A:,2A:+1^2A:+l,2A:+l] = (^2) 

The = {x^,x^)fi are the skew-scalar products of the monomials. By partial inte- 
gration, it is possible to derive a recursion relation relating t^j and 

Jo Jo 
= (/ + a - - ir (^i±it^) . (83) 

Using the antisymmetry of f^j = —fj^, all skew-scalar products can then be easily calcu- 
lated from ^0 = ''-i 1 = ^^id i • For a weight function with positive integer parameter 
a, the skew-scalar products are related to the case a = by t^j = t'^_^_^ j^^. Because of the 
antisymmetry of the skew-scalar product, (i?", R°;) = 0, but this relation does not impose 
an additional condition on the coefficients. 

We construct the skew-orthogonal polynomials from the homogeneous equations (|8TD . 
They can be easily normalized later by multiplying the even or odd order polynomials 
by a suitable scale factor. For the coefficients of the even order polynomials, the number 
of equations is one less than the number of coefficients, whereas for the odd order ones 
we lack two equations. To determine the polynomials, we fix 2k — ''^■zk+i 2k+i — 1 
and rf^_,_]^ 2fc = (^^e latter condition can be imposed because R2k+i is determined only 
up to a multiple of R2k)- In this way the polynomials can be determined accurately to 
about order 30. The skew-orthogonal polynomials for the Gaussian case can be derived 
in a similar way. 

To illustrate the asymptotic behavior of the skew-orthogonal polynomials we show in 
Fig. 1 the M-dependence of the ratio /fc(M)//jt(oo) (full curves) for a Gaussian (left) and 
a quartic (right) probability potential both with parameter a = . The integral Ik{M) is 
defined by 

/■(Af/^p(0))2 

4(M) = / dxe-'^'^^'^^Rkix). (84) 

Jo 
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Figure 2: The asymptotic behavior of the coefficients T2k,2k-i versus k for a quartic 
probabihty potential Its asymptotic value of —2 + -\/3 is depicted by the dashed curve. 

The weight function is given by 0a (a^) = x/2 for the Gaussian case and by (pai.^) = /2 for 
the quartic case. We have redefined M in units of 7rp(0), with p(0) the the average spectral 
density near zero (notice that with our convention for the weight function p(0) has a 
nontrivial /c- dependence). This figure shows that in the intermediate domain the integral 
/fc(M) oscillates around /fc(oo)/2. For even /c, the integral appears to converge in the 
oscillatory region. For odd k we show results for monic polynomials with normalizations 
r2k+i,2k = (middle figures) and by R2k+i{0) = (lower figures). In the ffist case, the 
odd order polynomials behave similarly to the even order ones, whereas in the second case 
the behavior is quite different, because the leading order asymptotic contributions cancel. 
The oscillations in the lower figures are still exactly about /fc(oo)/2 which in this case is 
close to zero. The even order skew-orthogonal polynomials always have as many positive 
zeros as the order of the polynomial. The odd order polynomials in the normalization 
'^2k+i,2k = are not very different from the preceding even order polynomials (see middle 
figure). However, typically one of their zeros is located on the negative real axis. For the 
normalization i?2A;+i(0) = 0, the total number of zeros is equal to 2k + 1, with one zero at 
X = 0. Fortunately, as we have seen in the previous section, the integral (^) over the odd 
order polynomials can always be tuned to zero so that we do not have to worry about the 
ambiguity of the asymptotic properties of the odd order skew-orthogonal polynomials. 

In Figure 2, we show the /c-dependence of the coefficients T2k,2k-i as defined in equation 
(0). Also shown is the analytical result for the asymptotic value of —2 + -\/3- The 
convergence to the asymptotic result is better illustrated by extrapolating to second order 
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Figure 3: The average spectral density p{u/Tip{Q)) /■npiOi) versus u. For a quartic proba- 
bility potential we show results obtained from the 30 lowest order skew-orthogonal poly- 
nomials (full curve) and the large-n analytical result (dashed curve). The semicircular 
distribution (dotted curve) represents the large-n average spectral density for a Gaussian 
probability potential. All results are for a = 0. 



in l/Zc. The values of T2k,2k—i 

at k ^ 10, 20 and 30 extrapolate to -0.26787 to be 
compared to -2 + V3 ^ -0.26795. 

In Figure 3, we show the average spectral density calculated from the quartic skew- 
orthogonal polynomials (full curve) and the analytical result (dashed curve) given by 



7rp(0) TT V 4n / V 8n 

For comparison we also show the semicircular distribution obtained for a Gaussian po- 
tential (dotted curve). 

In Fig. 4 we show the microscopic spectral density calculated from the first 30 skew- 
orthogonal polynomials for a quartic potential and a — 0. In the same figure we also show 
the result for a Gaussian potential. Clearly, the microscopic spectral density converges to 
the asymptotic result for n — > oo. Both in Fig. 3 and Fig. 4, the average spectral density 
p(0) depends on n because of the normalization of our weight function. 
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Figure 4: Microscopic spectral density Ps{u) for a Gaussian and a quartic probability 
potential. 

8 Conclusions 

We have shown universality for the chiral Orthogonal Ensembles. Our proof is based on 
a relation between the kernels for (3=1 and (3 = 2 and the universality of the kernel 
for (3 = 2. In this article we have completed the proof outlined in ||lO| by deriving an 
asymptotic property of the skew-orthogonal polynomials which relates an integral over 
the region near the largest zero to an integral in the microscopic region. Universality now 
has been shown for all three chiral ensembles. 

An alternative method for ensembles with (3 = 1 and (3 = 4 was recently proposed in 
This method does not rely on the construction of the skew-orthogonal polynomials 
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at all, but it is our point of view is that both methods are equivalent. For example, 
the matrix elements of some of the operators in our construction are also required in the 
method proposed by Widom. The advantage of our method is best illustrated by the 
result of this article in which universality has been proved by means of an asymptotic 
relation of skew-orthogonal polynomials. It would be interesting to identify this relation 
within Widom's approach. 

Finally, we wish to mention an alternative way of looking at universality. For theories 
with broken chiral symmetry and a mass gap we have two types of modes, the soft modes 
and the hard modes. An effective partition function is obtained by integrating out the 
hard modes. If the Goldstone bosons corresponding to the spontaneous breaking of chiral 
symmetry are the only soft modes, there is no need to do this calculation. The effective 
partition function can be written down solely on the basis of the symmetries of the theory 
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and is thus the same for all partition functions with the same global symmetries. The 
equivalence of the effective theory for the Goldstone modes and chiral Random Matrix 
Theory has been demonstrated nonperturbatively for the chiral Unitary Ensemble 
|60| , 0, |32|. For the the other two values of the Dyson index this has only been shown 
perturbatively |1T3| . 

The mass of the relevant Goldstone modes in the generating function of the spectrum 
is proportional to the square root of the distance to A = 0. Universal behavior is thus 
immanent in the microscopic limit. However, proving universality is equivalent to showing 
that the theory has a mass gap and that the Goldstone modes are the only soft modes. 
In QCD this is equivalent to proving confinement. 
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